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We study systems of coupled spin-gapped and gapless Luttinger liquids. First, we establish the 
existence of a sliding Luttinger liquid phase for a system of weakly coupled parallel quantum wires, 
with and without disorder. It is shown that the coupling can stabilize a Luttinger liquid phase 
in the presence of disorder. We then extend our analysis to a system of crossed Luttinger liquids 
and establish the stability of a non-Fermi liquid state: the crossed sliding Luttinger liquid phase 
(CSLL). In this phase the system exhibits a finite-temperature, long- wavelength, isotropic electric 
conductivity that diverges as a power law in temperature T as T — + 0. This two-dimensional system 
has many properties of a true isotropic Luttinger liquid, though at zero temperature it becomes 
anisotropic. An extension of this model to a three-dimensional stack exhibits a much higher in- 
plane conductivity than the conductivity in a perpendicular direction. 
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I. INTRODUCTION 

For over two decades a central theme in the study of 
correlated electronic systems has been the drive to un- 
derstand and classify electronic states that do not con- 
form to Landau's Fermi-liquid theory. A clear example 
of such "non-Fermi liquid" physics occurs in one dimen- 
sion (ID) Ip, where arbitrarily weak interactions destroy 
the Fermi surface and invalidate the notion of indepen- 
dent quasiparticles at low energy. Away from charge- and 
spin-density-wave instabilities, the ID interacting elec- 
tron gas forms a Luttinger liquid. The discontinuity in 
occupation at the Fermi energy is replaced by a power- 
law singularity, and the low-lying excitations are bosonic 
collective modes in which spin and charge decouple. 

In this paper, we investigate how such non-Fermi liq- 
uid behavior could arise in two- and three-dimensional 
systems consisting of arrays of quantum wires or chains. 
Our analysis is largely motivated by the unusual normal- 
state properties of High Tc materials. Prominent among 
these are: 

1. Resistivity perpendicular to the Cu02 plane much 
larger than the in-plane resistivity |^ . 

2. Linear temperature (T) dependence of resistivity 
along the conducting planes, and divergence 
of the Hall angle ||]. 

3. Angle resolved photoemission (ARPES) data show- 
ing a pseudogap and absence of dispersion in the c 
direction 0|. 

4. Linear temperature dependence of the Nuclear 
Magnetic Resonance (NMR) relaxation 1/Ti ||]. 

Anderson has suggested that these unusual normal-state 
properties of the cuprates are the result of non-Fermi- 
liquid physics in two dimensions The study of non- 
Fermi liquids in dimensions greater than one has, how- 



ever, proven to be quite difficult. Since the Fermi liq- 
uid is stable for weak interactions, perturbative methods 
about this state fail to locate non- Fermi liquid states . 
Moreover, generalizations of the bosonization technique 
to isotropic systems in higher dimensions have indicated 
that Fermi-liquid theory survives provided the interac- 
tions are not pathologically long-ranged [|| . 

An alternative approach has been to study anisotropic 
systems consisting of arrays of parallel weakly coupled ID 
wires ^. Coupled Luttinger liquids have been studied 
extensively for the past thirty years, mostly in the con- 
text of quasi one-dimensional conductors. Single-particle 
hoppings as well as pair-hopping correlations were shown 
to destabilize the Luttinger liquid phase. Following An- 
derson's suggestions regarding high-temperature super- 
conductors, there have been a series of studies on cou- 
pled quantum chains, which confirm this view [ p^ though 
the issue remains controversial JlT| . It has recently been 
proposed p^ , p^ that for a range of interwire charge and 
current interactions, there is a smectic-metal (SM) phase 
in which Josephson, charge- and spin-density-wave, and 
single-particle couplings are irrelevant. This phase is 
a two or three-dimensional anisotropic sliding Luttinger 
liquid whose transport properties exhibit power-law sin- 
gularities like those of a ID Luttinger liquid. It is the 
quantum analog of the sliding phases of coupled classi- 
cal XY models first identified in the context of phases of 
DNA-lipid complexes 1 14 1^ . The study of coupled Lut- 
tinger liquids becomes particularly relevant in the con- 
text of striped phases which have recently been found 
in quantum Hall systems |l^ and in the cuprates [ pT[ . 
These phases are characterized by inhomogenious dis- 
tributions of charge and spin, where, for example, the 
charge carriers might be confined to separate linear re- 
gions, thus resembling stripes. Sliding Luttinger Liquid 
phases might also be relevant for a variety of other sys- 
tems such as quasi one-dimensional organic conductors 
and ropes of nanotubes. 
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The work on sliding Luttinger liquid phases has also 
been extended to a square network of ID wires formed 
by coupling two perpendicular smectic metals , which 
has been shown to exhibit a crossed sliding Luttinger Liq- 
uid (CSLL) phase. At finite temperature T, the CSLL 
phase is a 2D Luttinger liquid with an isotropic long- 
wavelength conductivity that diverges as a power-law in 
T as T ^ 0. At T = 0, it is essentially two independent 
smectic metals. This model could be realized in man- 
made structures constructed from quantum wires such 
as carbon nanotubes. Extension of the model to a three- 
dimensional stack may be relevant to the stripe phases 
of the cuprates. Based on neutron and x-ray scatter- 
ing measurements, it has been suggested that spin-charge 
stripes in the adjacent Cu02 plane are orthogonal to each 
other 0. 

In order to understand the nature of the sliding phase, 
it is easiest to think in terms of the classical analog. 
Imagine a stack of 2D XY models coupled by the Joseph- 
son coupling cos{9n—0n+i) where 9n is the XY-angle vari- 
able in layer n. Such a system always goes directly from 
a 3D ordered phase to a completely disordered phase as 
a function of the temperature T. However in the pres- 
ence of interlayer gradient coupling terms of the form 
V^m • V^n, the system may have an intermediate slid- 
ing phase [|l5| , that exhibits an in-plane 2D order with 
power-law decay of correlations along the planes. In the 
absence of Josephson couplings, in the sliding phase, 9 
in neighboring layers freely slide over each other with 
no energy cost, and two-point correlation functions are 
identical in form to those of a stack of decoupled 2D lay- 
ers. When inter-layer Josephson couplings are present, 
though irrelevant, two-point correlation functions decay 
exponentially perpendicular to the layers. Thus the slid- 
ing phase has in-plane 2D order but is disordered in the 
third direction. For sliding Luttinger liquid models, the 
2D layers are replaced by ID quantum wires, and inter- 
wire current-current and density-density couplings play 
the role of the gradient couplings. Thus, in the sliding 
Luttinger liquid phases, the correlation functions along 
a given wire exhibit the same power-law functional form 
as in a ID Luttinger liquid. 

In this paper, we review and discuss in detail results on 
the sliding Luttinger liquid phase. Our discussion of the 
parallel Luttinger liquids follows Ref. jl^], but it estab- 
lishes the stability of the sliding phase to a more complete 
set of inter-wire operators and also to disorder. For the 
crossed sliding Luttinger liquid phase (CSLL), we provide 
detailed calculations that were presented only briefly in 
Ref. jl^. The stability and transport properties of the 
CSLL phase form the central results of the paper. We 
also extend the analysis to a three-dimensional stack of 
crossed Luttinger liquids, and show that we can still ob- 
tain a sliding Luttinger liquid phase, with hopping be- 
tween planes being irrelevant. In this phase, the in-plane 
and perpendicular conductivity can be quite different and 
exhibit different temperature dependences. This could be 
of relevance to the normal state of the cuprates. 



This paper is organized as follows. In Sec. II we ex- 
plore the stability of the sliding phase for an array of 
parallel coupled Luttinger liquids with respect to a large 
(but not complete) set of inter- wire operators. For a one- 
dimensional system of interacting spin-1/2 fermions, the 
spin excitations could either be gapped or gapless. In 
the spin-gapped Luther-Emery regime, the system can 
by described by a single Luttinger liquid for charge. In 
the gapless case, both spin and charge are dynamical de- 
grees of freedom, and there are two Luttinger parameters 
(kct, tp), and two velocities {va-, Vp). We study coupled 
Luttinger liquids both for the spin-gapped and the gap- 
less case, and in each case demonstrate the stability of 
the sliding phase. In addition, we study the effect of dis- 
order. We find that density-density and current-current 
interactions that stabilize the sliding phase, also make 
disorder more strongly irrelevant. Note, however, that 
in this paper we do not establish stability with respect 
to all multi-wire operators. Since the overall strength of 
these higher order interactions that we neglect are ex- 
pected to be much smaller than those we consider here, 
their relevance becomes important only at very small 
temperatures. We delay a more complete study of their 
effects to a future publication. Throughout the paper 
we will use the term "stable" somewhat loosely to refer 
to stability with respect to a restricted set of interac- 
tions that include all two-wire interactions. In Sec. Ill 
we extend the above analysis to two identical coupled 
Luttinger liquid arrays, arranged such that the wires of 
one array run perpendicular to those of the other. In 
this case we establish the stability of a crossed sliding 
Luttinger liquid (CSLL) phase to a large class of opera- 
tors. The inter-array density-density couplings effectively 
renormalizc the intra-array couplings; however, the sta- 
bility of the CSLL phase turns out to be identical to 
the stability of the smectic-metal phase for each array, 
but with the intra-array couplings replaced by the renor- 
malized couplings. We then generalize our analysis to a 
three-dimensional stack of arrays, with wires on each ar- 
ray running perpendicular to the wires of the consecutive 
array, and obtain once more a sliding phase that appears 
stable. In Sec IV we explore the transport properties 
of sliding phases. In particular, we focus on the power- 
law singularities of the conductivity as a function of T 
as T — > 0. Finally Sec V sums up our principal results. 
Appendix A sketches out the steps used to obtain intra- 
wire correlation functions, Appendix B and Appendix C 
present details of some integrals used in Sees. II and III, 
and Appendix D presents details of conductivity calcula- 
tions using the Kubo formula. 

II. COUPLED PARALLEL LUTTINGER LIQUIDS 

Our goal is to construct non-Fermi-liquid electron sys- 
tems in two and three dimensions. Our basic building 
blocks are one-dimensional quantum wires that exhibit 
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Luttinger-liquid phases with non-Fermi hquid power-law 
decay of correlation functions. We couple these wires to- 
gether in arrays in such a way that they retain their one- 
dimensional non-Fermi liquid character yet allow non- 
vanishing intcrwirc electron transport at nonzero fre- 
quency or temperature. We will consider several types of 
arrays, and it is useful to introduce a nomenclature for 
them. The long axis of parallel wires defines one direction 
in space. The wires can be arranged in either a one- or 
a two-dimensional Bravais lattice of points in the plane 
perpendicular to that direction. We, therefore introduce 
the notation d±:l to denote an array of parallel wires cen- 
tered on lattice points in a dj^-dimensional space. The 
resulting structure occupies a dj^ + 1 dimensional space. 
Thus a 1:1 array is a planar array of equally spaced paral- 
lel wires, and a 2:1 array is a three-dimensional columnar 
array of wires. As stated, we will restrict our attention to 
arrays in which the wires lie on a dj^-dimensional periodic 
lattice. For 2:1 arrays, we will only consider, in detail, 
those based on a 2-dimensional rectangular lattice. Both 
1:1 and 2:1 arrays can exhibit anisotropic sliding phases 
for an appropriate choice of interwire potentials. 

In this section we will consider anisotropic sliding 
metal phases in 1:1 and 2:1 arrays. In the succeeding sec- 
tions, we will discuss how crossed 1:1 and 2:1 arrays can 
produce two and three-dimensional sliding phases with 
Civ symmetry and isotropic long wavelength finite tem- 
perature conductivity. We begin our discussion with a 
review of ID Luttinger liquids. 



A. Review of ID Luttinger liquids 



)/V2, 



/2, (2.3) 



and similarly for the (jj variables. Here p characterizes 
the charge excitations, and a the spin excitations. For 
a single Luttinger liquid, we can write down an effective 
Hamiltonian for the charge sector 



(2.4) 



where Vp is the velocity, and Luttinger liquid pa- 

rameter (this is the inverse of the usual Luttinger liquid 
parameter g). Alternately, we could write down the ac- 
tion, in a path-integral formulation, as 



-2l{d.Jp){dr(l^p)]. 



Integrating out either the Op or the 
tain 



K.p{dx(l)pf 

(2.5) 

<bn variables, we ob- 
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In the spin sector, we might either have a spin gap cor- 
responding to the Luther-Emery regime, or have a gap- 
less phase where the spin excitations are described by 
a Hamiltonian of the same form as Eq. (2.4) with the 
parameters and Va- In the gapless phase, SU(2) sym- 
metry imposes the constraint k.^ = 1. 



One dimensional Luttinger liquids are most easily de- 
scribed in terms of collective bosonic charge and spin- 
density modes. The bosonized form of the fermion oper- 
ators near the left and right Fermi points are p9| 




(2.1) 



where R and L stand for the right and left moving elec- 
trons, respectively, s is the spin index, j denotes the wire 
number, e is some intra-chain cutoff, and Tjn/j^ g j are the 
Klein factors. We can write down an effective theory for 
the low-energy excitations in terms of the boson opera- 
tors $. It is convenient to define 



($fl,,sj + «>L.sj)/V4^, 
(*ii,sj - $L,sj)/%/4^, 



(2.2) 



where 9s j is a phase variable, and (psj is the conjugate 
density variable pO| . Since the bosonic excitations can 
be characterized in terms of charge and spin excitations, 
we further define 



B. The spin-gapped 1:1 array 

In this subsection we consider the simplest array of 
quantum wires, the two-dimensional 1:1 array of Lut- 
tinger liquids in the spin-gapped phase. It has been sug- 
gested that this case might describe the striped phases of 
high-temperature superconductors pl| . In this subsec- 
tion all bosonic variables refer to the charge sector, and 
we do not write explicitly the subscript p. In general, we 
expect density-density and current-current interactions 
between the wires, which can be represented by a action 
of the form 

•Sint = ^ J dxdTjp_^nix,T)Wf_i,{n-n')jf_,^ri'ix,T), 



n,'n' 



(2.7) 



where j^,„ = (p„(a::, r), J„(a;, r)) with p„ = dx(j)^n{x,T) 
the density and J„ = OxOnix^r) the current on the nth 
wire. The density-density interaction is an effective in- 
teraction generated by both the screened Coulomb and 
the electron-phonon interaction. In the striped phases, 
stripe fluctuations lead to current-current interactions 
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pl[ . These current-current and density-density interac- 
tions are marginal and should be included in the fixed- 
point action. They are invariant under the "sliding" 
transformations — > + a„ and 0„ — > 0„ -I- . Equa- 



tions (2.4) and (2.7) define the fixed-point action of the 
smectic- metal phase [|2| , which can be written as 



n+j 



J2 Vfid-Mid.Cpn+j) + 2i{d^9n)idr<l>n)]. 



(2.8) 



Upon integration over or On, respectively, the effective 
action for 0„ and d)„ become 



So 



1 



(PQ 1 
(2^2 
d^Q 1 



1 



v{q^)q^\\em' (2.9) 



where Q = {uj, q\\, qj_), with g|| the momentum along the 
chain and q± that perpendicular to the chains. Here 



n{q±) = ^VHq±)/V%q^) 
= yJvHq±)VO{q±), 



(2.10) 



p = = + L^L + [L^R + c.c] 



(2.11) 



Thus the density has two pieces: the first piece is the 
course-grained density and the second is modulated at 
2kp where kp is the Fermi wave-vector. Correspond- 
ingly, the density-density correlation has two pieces. We 
consider G^{x,t): the component of the density-density 
correlation function with 2kF modulation. Thus 

G^{x,t) = (e'V27('^'.(=^^^)-<^.(0'0))+'2/cr^) +C.C.. 
It easy to show that for large x (see Appendix A), 
Ai cos{2kpx) 



where 



G^ix,0) 



AcDW,cc- = 



dqi 



1 



^ 2tt Kiqj_) 



(2.12) 



(2.13) 



FIG. 1. A schematic depiction of a 2-dimensional array of 
quantum wires 



and Ai is a constant. Alternately, for large r 



(2.14) 



where A2 is a different constant. In general, Gtj){x, r) can 
be written in the scaling form 



G^{x,t) « X 



-Ac 



(f) 



where V^{qA_) and V'^{qA_) are the Fourier transforms of where 

and Vj" with respect to the wire index j. Eqs. (2.8) 
and (2.9) define the action of the ideal 2D sliding Lut- 
tinger liquid or SM phase. Even though they include 
interwire interactions described by Eq. (2.7), they yield 
power-law correlations characteristic of a ID Luttinger 
liquid. 

We could consider, for example, correlation functions 
involving the density variable p. We note that 



F{y) Ai as y ^ 



as y — > 00. 



(2.15) 



(2.16) 



A variety of interactions, other than those of Eq. (2.7) 
couple neighboring wires. The sliding phase is stable only 
if the interactions are irrelevant, that is, only if they scale 
to zero with system size. We will now investigate pertur- 
batively the relevance of these interchain interactions to 
determine under what conditions the sliding phase is sta- 
ble. Due to the spin gap, single-particle hopping between 
chains is irrelevant, and the inter-chain interactions that 
could become relevant are the Josephson (SC) and CDW 
couplings, which are represented by operators of the form 



X(i?tj+„ij.j+n + i?J.j+„L|j+„)] + C.C, (2.17) 



and 



x(4j+n^Tj+« + Llj+n^hj+n)] + CC. (2.18) 

At low temperatures the spin variable (f)^ is effectively 
frozen, and these interactions depend only on 9i. Their 
associated actions can be expressed as: 



4 



5'sc,n ^ J dxdr cos[\/27r(6'i - 9i+n)\, 



ScDW,n = V,i 



dxdr cos[v 27r( 



(/.,+„)]. (2.19) 



The relevance of these terms are determined by the 
scahng dimensions of the corresponding operators, 
cos[\/2n{di — Oi+n)] and cos[\/27r(0i — which are, 

respectively. 



Asc,n = J "^(1 ~ cos nq_i_)K{q_i_) 
^ dq± (1 — cosnq_L) 



AcDW.n / r, f \ 

' 2tt K[qj_) 
The exponent Agc,ri follows from 



(cos v27r( 



)) 



(2.20) 



(2.21) 



and 
1 



d^'Q 



(27r)3 
" dq± 



(\eiq,u)\^){l~cosiq^n)) 

{k((7j_)(1 - cosg_Ln) 
dquduj 1 



}. (2.22) 



where () denotes averaging with respect to S in Eq. (2.8). 
The integral in the square brackets diverges logarithmi- 
cally with system size L (~ ClnL) Using this, we 
find 



(2.23) 



(cos V27r(6'j - 9j+n)) ^ i" 



where Asc.n is given by Eq. (2.20). From the above 
equation, it follows that 



SC,n) 



(2.24) 



Similar calculations produce Acdw.oo- For a stable smec- 
tic metal phase, these terms have to be irrelevant, imply- 
ing 



A 



CDW,n 



> 2, 



> 2 



(2.25) 



for all n and n . 

If any Agc,„ < 2, the smectic metal (SM) phase is 
unstable to the formation of an anisotropic 2D supercon- 
ductor. If any AcDW.n < 2, the SM phase will flow to 
a 2D longitudinal CDW-crystalline phase with 2k p den- 
sity modulations along the wires and phase locked from 
wire to wire. Notice that if K{q±) is uniformly small, 
Asc.n's are small and AcDW,n's are large, and for large 
K, Asc.n's are large and AcDW,n's small. For a stable 
sliding phase we need all Asc,ra's and AcDW,n's to be 
greater than two. Our strategy to create a stable sliding 




FIG. 2. K{q) as a function of q/n, with a minimum, Kmin 
at g = (Jo- 



phase is to choose K{q±) of the form shown in Fig. (2), 
with K having a minimum, Kminj at q± = qo- When Kmin 
becomes zero, the system undergoes a transition to a 
transverse CDW modulation with wavevector qo, where 
the charge-density varies from wire to wire. For small 
but positive Kmin, the system is close to the transverse 
CDW instability. As pointed out in ||l^, a transverse 
CDW would frustrate crystallization of fcrmions, since 
kp is now a function of the chain index j. Strong fluctu- 
ations of this kind prevent the locking in of density fluc- 
tuations along the wires. Thus, in order to stabilize the 
sliding Luttinger liquid phase, we need to tune Kmin to be 
very small compared to average K{q±), so that Asc,n and 
AcDW.n can both be made large. Note that in addition 
we should also consider interchain operators of the form 
i?! ^i?! ■R^^jj^nRi,j+n- These interactions, however, turn 
out to be automatically irrelevant if the superconducting 
and CDW interactions are irrelevant, and hence merit no 
further consideration. 

In addition to the lowest-order interactions between 
pairs of chains described by Sa,c.n and S'cDW.n, there 
can in principle be higher-order multi-chain interactions 
with actions of the form 

Sscsp = X! / dxdTjs^cos[\/2Ti(^Spe.i+p)], 

•ScDW^Sp = X! / ^2;dTV,5j, cos[V27r(^ Sp0j+p)] (2.26) 
i ■' p 

where jTsp are the inter-chain Josephson couplings, V^j, 
the inter-chain particle-hole (CDW) interactions, and 
Sp is an integer-valued function of the chain number 
p satisfying Sp = 0. The scaling dimensions of 

cos[V2n(J2p SpOi+p)] and cos[V2n(J2p Spipi+p)] are, re- 
spectively. 
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iCDW,s„ 



27r 



^SpSp' cos{{p~p )q±} 
p,p' 

(2.27) 



These perturbations are irrelevant if 

AcDW,sp > 2, Ag(^ > 2 



(2.28) 



for all sets of Sp and 



The relevance of higher order 
terms of this form is a subtle issue, and we will elaborate 
on this in a future publication. However the strength of 
these terms, as measured for example by , is small and 
they become important only at very small temperatures 
even if they are relevant. In this paper we will restrict 
ourselves to pairwise CDW couplings of the form given 
by Eq. (2.17); however, we comment on the higher order 
superconducting terms below. 

To explore the regions of stability of the smectic metal 
(SM) phase, we follow Refs. and (l5| and take 

K{qi_) = K[l + Ai cos{qi_) + Aa cos(2q_L)]- (2.29) 

The parameters Ai and A2 can be tuned to set the value 
9o of g_L at which k((7j_) reaches a minimum and the min- 
imum value of K.{q±_), K{qo) = KA. For a specific value 
of A and k, 



Ai = 
A2 = 



4(1 - A) cosgo 
~(l + 2cos2((7o)) 

(1-A) 
(l + 2cos2go)' 



(2.30) 



unless go = or tt, in which case only A1+A2 is fixed by A. 
We note that in the above equation (i.e. for qg ^ 0, tt), 
while A2 is always positive for A < 1, Ai can be either 
positive or negative. Typically positive A's correspond to 
repulsive interactions. 

In this paper we treat K as the control parameter. For 
small K the system becomes superconducting while for 
large K the system goes into a CDW crystalline phase. 
We will show that through judicious tuning of Ai and A2 
it is possible to have an intermediate window of K where 
the smectic metal phase is stable. For this purpose, we 
define a^^ — Asc,s„/-?^ and 6„ = AcDW,n-^, where 
and bn depend only on Ai and A2, and 



^^^(4 + «p(''^p+i + Sp-i)Ai/2 



+ Sp{Sp+2 + S)9-2)A2/2). 



(2.31) 



The SM phase becomes unstable to inter-chain Joseph- 
son couplings for K less than K^c = maxs^ (2/asj,) and 
unstable to inter-chain CDW interactions for K greater 



than i^cDW = iniiin(&n/2). Thus the smectic metal 
phase is stable over a window of K, Kgc < K < Kqbw, 
provided 



(3: 



Ksc 



> 1. 



(2.32) 



min.wrt.m&Srj 



We note, once more, that in this paper we consider sta- 
bility with respect to all superconducting terms, but only 
pairwise CDW terms of the form given by (2.17). If 
< 1, the system goes directly from a 2D supercon- 
ducting (SC) phase to a CDW crystal as K is increased, 
without passing through the smectic metal (SM) phase. 
For a stable sliding phase, we need to make A very small. 
The value of AcDW.n for A small is determined by values 
of q± near go- We can therefore set K{q±_) 



li{q^) « K[A + C{qi_ - qnf 



where 



^^^^2A2sin2go. 



This gives us 



A 



CDW.n 



2K 



K{\ — cos(nqo)e 



(2.33) 



(2.34) 



(2.35) 



where C has been defined in Eq. ( 2.34 ). We set A = 
10^^. We consider the range of qo/T^ lying between 0.25 
and 0.75. This range can be broken into three sections: 

1. 0.25 < go/TT < 0.41957. 

In this range Ai < and |Ai| > A2. Here the 
most relevant superconducting term corresponds to 
the multi-chain operator cos[-\/27r(6'i + 9i+i — Oi+z — 
Oi+4)] 1^^. The dimension of this operators sets 
the minimum of . Thus, in this range, minos = 
(2 + Ai - A2/2) and Ksc = 1/(2 + Ai - A2/2). ' 

2. 0.41957 < qo/n < 0.5804. 

In this region |Ai| < A2. We find that is smallest 
for the set s„ ~ Sn q — Sn 2- Thus, in this range 
T^sc = 2/(1 - A2/2); 

3. 0.5804 < qo/n < 0.75. 

Here Ai > A2, and as is the smallest for the set 

Sn = Snfi - 5n,l. ThusVgc = 2/(1 - Ai/2). 

In Fig 3 we plot /3 as a function of go- The minima of 
the curve correspond to go — 27rZ/m, where I and m are 
integers. Also note that since Ai has the same sign as 
(—cos go), there are regions of stable smectic phase for 
positive as well as negative values of Ai. 

Having established a stable smectic phase for the pure 
system, we now study the relevance of quenched disorder 
in this phase. Disorder gives rise to a random electron 
potential, D{x), with associated action 
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(2.36) 



D(x) can be treated as a Gaussian random variable, with 
zero mean and local fluctuations such that 



D{x) = 



D{x)D{x') = 1^d6{x -x). 



(2.37) 



where the over-line signifies averages over the random- 
ness. By a generalization of the Harris criterion |Q, it 
can be shown quite easily that (also, see Giamarchi and 
Schulz [p5|) disorder is irrelevant if Acdw,oo > 3, where 



A, 



CDW.oo 



dq 1 

-TT 27r K.{q) 



1 1 



(2.38) 



For the range of parameters we are considering where the 
SM phase is stable, Acdw,oo is large and thus disorder 
is strongly irrelevant. This is an important point. For a 
single Luttinger liquid in the repulsive region, k > 1 and 
disorder is always relevant. However inter- wire interac- 
tions can drive disorder irrelevant for k((?_l =0) > 1, 
even in regions of phase space where all interactions are 
repulsive. 

Thus, there is a small but finite region of phase space 
where the smectic metal phase appears stable. We should 
note that over a larger region of phase space the only rel- 
evant operators involve nonlocal interactions of the form 
V„ cos[\/27r(0i — (pi+n)], where V„ is expected to be ex- 
ponentially small, for large n, in the bare Hamiltonian. 
Though relevant, these operators would only play a role 
for fcsT smaller than some energy scale set by V„. So, 
for example, there will be a range of temperatures where 
we will only need to consider the relevance of Vi and Ji . 
These can be made irrelevant over a reasonably large re- 
gion of phase space (see [|l2|). Thus, even though the 
region of phase space where the smectic phase is strictly 
stable is highly restricted, at finite temperature and for 
weak coupling, we expect a much larger region of phase 
space whose behavior is governed by the sliding Luttinger 
liquid ground state. 



C. The gapless 1:1 array 

We now consider 1:1 arrays of wires in which both 
charge and spin excitations are gapless. In this case, 
there are two Luttinger liquid parameters (Kp, k^) for 
the charge and spin modes respectively, and two veloci- 
ties {vp,Va) on each wire. To maintain gapless Luttinger 
liquids and SU(2) spin symmetry, we do not include any 
marginal spin-spin coupling terms in the Hamiltonian. 
Thus the spin degrees of freedom are represented by the 
fixed-point action 
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(2.39) 




0.25 



0.5 



0.75 



FIG. 3. Plot of /3 = KcDw / Ksc as a function of of go/vr. 
For /3 > 1 there exists a region of K over which the non-Fermi 
liquid phase is stable. 

with Kcr = 1. In a more general treatment, one could 
include spin-spin coupling terms and consider their rele- 
vance/irrelevance, maintaining, however, the SU(2) sym- 
metry of the spin-sector. We leave that for future con- 
sideration. The charge modes are still represented by 
Eq. (2.9), with K{q±) and v{q±) replaced by Kp{q ±) and 
Vp{q±). The form of Kp is still given by Eq. ( p.29[ ). 

We again consider the relevance of single-particle, 
CDW and SC tunnelling. The SC and CDW tunnel- 
ing were already considered in the previous subsection. 
When the spin variables are included, Eqs. (2.11) and 
(2.12) become 



SsG,n ^ Jn dxdr y^ cos[ 

3 



/27r( 



n)] 



ScDW,n 



cos( 



2tt 

dxdr 



crj) COs(V27r0CTj+n) 



3 



cos[v27r( 



^P,3+' 



cos(v 27r0crj) cos(v27r(?!)CTj+„) 



(2.40) 



The 0O- variables now contribute to the dimensions of 
these terms. Because k„ is constrained to be 1, the con- 
tribution of the a variables is trivial, and the dimensions 
of these terms are given by 



, (gap) 



1 



AcDW,n — A 



(gap) 
CDW.n 



(2.41) 



where A||p^ „ and Ag^'^^ are given by Eqs . (2.13) and 
(2.14) with K replaced by Hp. 

Since the a variables are no longer gapped, single- 
electron tunneling is no longer irrelevant. Single-particle 
hopping is described by operators such as ^j,t^]+„|j 
which can be represented by terms of the form 



, (gap) 
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X < e 



{■ 



(2.42) 



The expectation value of the term in the curly bracket 
goes as L~^/^ as system size L goes to infinity. Thus 



^SC,nJ 



(2.43) 



Regions of phase space where Ssc.n is relevant corre- 
sponds to the superconducting phase, whereas regions 
of relevance of ScDW,n correspond to the CDW crys- 
tal phase. Regions where both these are irrelevant, 
but single-particle hopping is relevant correspond to the 
Fermi-metal phase. For a stable smectic metal phase, we 
require that all these operators be irrelevant. The su- 
perconducting and CDW coupling terms are irrelevant 
if 



. (gap) 
CDW.n' 



> 1, 



(2.44) 



for all n and n The condition for single particle hopping 
to be irrelevant is that 



A (gap) r. 



(2.45) 



for all n. 

We now proceed exactly as for the gappe d cas e, as- 
suming Kp to have a form as given by Eq. ( 2.29| ). As 

before, we may write A^f^^ , 



n IK a'*^''^^ ~ h K 



Ksc — niax„(l/a„) and i^cDW = min„(6„). Provided 
KcT>w / K'=ic > 1, there is a window of K, Ksc < K < 
KcDW where the system is stable with respect to both 
the CDW and superconducting couplings. For the smec- 
tic metal phase to be stable, the single particle hopping 
has to be irrelevant as well, which indicates that 



K 



b„K > 6. 



(2.46) 



This condition is violated for K lying between and 
_ftr_|_ where = min„_ftr_^„ and = max„iir+^„, with 



K. 



3 ± V9 - anbn 



±,n 



(2.47) 



The single electron hopping is relevant in a large region 
of phase space, indicating an instability towards a Fermi 
liquid (FL) phase . We write Ai and A2 as functions of A 
and qq, (see Eqn. 2.30), and set A = 10~^. Higher order 
terms involving 6^ and (jj^ in general are less relevant in 
this case, and we do not need to consider the whole set of 
operators j2^. Depending on gg, we have the following 
possibilities for phases as K is increased : 

1. Superconductor(SC) -> Fermi hquid(FL)^ CDW 
crystal; 



0.69 



o 

0.67 



0.65 





SM 


sc 
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FL 


CDW 
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FIG. 4. A plot of the phase diagram in go, -ft" space with 
A = 10^^. SC stands for Superconducting, FL for Fermi 
liquid, SM for smectic metal, and CDW for charge density 
wave crystal. 



2. SC ^ Smectic Metal(SM) CDW crystal; 

3. SC ^ FL ^ SM ^ CDW crystal; 

4. SC SM ^ FL -> SM ^ CDW crystal. 

The phase diagram is complicated, and we plot a re- 
gion of K, go space in figure 4, in the absence of dis- 
order. Backscattering due to disorder is irrelevant for 
AcDW 00 > 2' which is automatically satisfied in the SM 



phase. 



D. The three-dimensional anisotropic sliding phase 

We now turn to three-dimensional 2:1 arrays with wires 
on a periodic 2D lattice with primitive translation vec- 
tors fli and (12 ■ Each wire occupies a position riioi -I- 71,2(12 
on a 2D lattice and is labeled by the integer valued vec- 
tor n = (m, 712). We will focus on the spin-gapped case, 
though extensions to the gapless case proceed exactly as 
in the previous subsection. Taking into account inter- 
wire density-density and current-current interactions we 
could write down an action of the form (2.7), but where 
we now have sums over columns in a 2D lattice, with 



jfj.,n{x,T) and n — {rix^ny). When trans- 



formed to Fourier space, this action becomes 



S 



[V'{qi^)qt\OY 



(2^)4' 

+V^{q^)q}M\'' -iLoq\\{e*cj>., 



+ C.C.} 



(2.48) 



where Q = (a;,g||,g_L) with q±^ a vector in the first Bril- 
louin zone of the 2D lattice of columns. We choose the 
X-axis to lie along wires, so that q±^ = (qy,qz). The 9 
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FIG. 5. A 3-dimensional array of quantum wires. 



or the <j) variables may be integrated out, giving us the 
effective actions 



Se 



Sik = 



where 



(27r) 



4 9«(9-L) 



1 



(2^)2 2Ac(ql) 



1 



+ t;(g-l)9f||^(Q)|^ (2.49) 



v{q±) = ^Jvf{q,_)V6{q^). 



(2.50) 



In three dimensions it turns out that stabiUty of the shd- 
ing phase with respect to the complete set of operators 
requires an even further fine-tuning of the generalized 
current-current coupling terms. In particular, k(%,(?z) 
should have a minimum KA at some qy = qo^y, qz ~ qo,z, 
with both A and the second derivative of k/K being 
much smaller than unity at the minimum. Let us con- 
sider two examples of the form that K{qy,qz) could as- 
sume in order to obtain a stable sliding phase. 

The first example is one which is symmetric with re- 
spect to qy and g^. We assume that the wires are ar- 
ranged in a square or rectangular pattern, and align the 
y- and z- axes along the edges of the rectangle. We con- 
sider tha form: 

«^(9y. qz) = K[l + Ai cos{qy) + Ai cos{qz) 

+X2COs{qy)cos{qz)]'^. (2-51) 

Ai and A2 are adjusted such that k has a minimum KA^ 
at qy = qz = qo- This gives 



Ai = - 

A2 = 



cos (go) ' 
(1-A) 



cos2(go)' 

Close to {qy, qz) = {qo, qo), we can expand k as 

K ~ K[A + X2{qy - qo)(qz - qo)]"^- 



(2.52) 



(2.53) 



As before Asc,m = Kam and AcDW.n = Kbn, where 
m and n are now vectors. The sliding phase is stable 
provided 



Ksc ~ 



> 1. 



(2.54) 



min.wrt.m&n 



For A = 10^2 (or smaller) there is a large range of qo 
where the sliding phase is stable. 

One could also consider the highly anisotropic form 



K{qy, qz) = K[{1 + Ai cos{qy) + A2 cos{2qy))'^ 
+X3{1 + cos{qz))]. 



(2.55) 



Again, for any A3, one can adjust Ai and A2 to produce 
a stable sliding phase. We conclude by noting that in 
three dimensions obtaining a sliding phase requires an 
even finer adjustment of parameters than in 2D. At fi- 
nite temperature, as before, the region of phase space 
controlled by the smectic metal fixed point is expected 
to expand considerably. 



III. CROSSED SLIDING LUTTINGER LIQUID 
PHASE 

Having established regions of stability of the sliding 

metal phases formed from arrays of quantum wires, we 
now turn to the investigation of sliding phases formed 
from crossed arrays of wires. We consider two basic con- 
figurations: one a two-dimensional system formed from 
two coupled 2D sliding phases (1:1 arrays) oriented at 
right angles to each other and the other a 3D system 
formed by stacking the crossed two-dimensional system. 
The latter three-dimensional system can be constructed 
from two inter-penetrating 3D anisotropic sliding phases 
(2:1 arrays), of the type discussed in Sec 2.D, oriented at 
right angles to each other. Both the 2D and 3D sys- 
tems have symmetry. As a result, their in-plane 
conductivities at finite temperature are isotropic at long- 
wavelengths. In this section we demonstrate the exis- 
tence of a sliding phase in the crossed arrays that is sta- 
ble if the sliding phase in the constituent arrays is stable. 
The correlation functions in this phase exhibit power-law 
decay along the planes, and the electric conductivity di- 
verges as a power-law in temperature T as T — > 0. 
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FIG. 6. Schematic depiction of a 2-dimensional crossed ar- 
ray. 

A. Crossed 2D sliding phase 

We consider now a square grid of wires, starting again 
with the spin-gapped case. The system consists of two 
arrays of quantum wires, the X- and Y- arrays running, 
respectively, parallel to the x- and y-directions. Each 
wire sees a periodic one-electron potential from the ar- 
ray of wires crossing it. For simplicity we assume that 
this periodicity is commensurate with bands in the wire. 
This leads to a new band structure with new band gaps. 
It is assumed that the Fermi surface is between gaps so 
that the wires would be conductors in the absence of fur- 
ther interactions. By removing degrees of freedom with 
wavelengths smaller than the inverse wire separation, we 
obtain a new effective theory whose form is identical to 
the theory before the periodic potential was introduced. 
Thus, in the absence of two-particle interactions between 
crossed arrays, the system could be in a phase consisting 
of two crossed, non-interacting smectic-metal states. 

We will now demonstrate the existence of stable sliding 
phase in the crossed arrays. In addition to the interwire 
couplings within each array, we need to consider Coulomb 
interactions between wires on the X-array and wires on 
the Y- array. Those inter- array couplings arc marginal 
and should be included in the fixed point. They do not, 
however, change the dimensions of the operators, except 
by rcnormalizing K(q±). For a stable sliding phase, ad- 
ditional interactions between the two arrays, such as the 
Josephson and CDW couplings, have to be irrelevant. We 
will show that it is possible to tune K{q±) such that this 
is indeed the case. 

The Coulomb interactions between electrons on inter- 
secting wires give rise to a term in the Hamiltonian of the 
form y^ „(a;,?/)/92:,m(a;)/9y,„(y), where (a;) [py,m{y)] is 
the electron density on the mth wire on the X{Y)-aiTay 
at position x {y). We expect ^{x, y) to have the form 
V'^{x — na, y — ma), where a is the distance between par- 
allel wires. Thus we represent the interaction between 
the X and F-array as 



/ 



dxdy[d:c(j)x,m.V{x -na,y- mb)dy(j 



If all parameters for the X and Y-arrays are the same, 
the crossed-grid action as a functional of the 9 and 
variables can be written as 

+ViQx, qy)<lxqv{(t>x4'*y + c.c.} 

-iujqx{9l(t>x + c.c.} - iujqy{ei4>y + c.c.}] (3.1) 

with obvious definitions for (p^ = 4'x{'-^,QxjQy)y4'y^(^xi 
and 6y. It should be noted that this is an effective theory 
with — - < qx,qy < —■ Integrating out the (f) variables, 
we are left with an effective action, which is conveniently 
expressed in matrix form as 

Se = \j d^kdwea{G-%j,0t (3.2) 

where a = x,y; b = x,y. Here 



G 



(3.3) 



where 



Kx(q) = 
^^x(q) = 



7(q) 



VHqxW^iqy)' 
V^g,)7(q) 



V V<^{qx) 



7(q) 

7(q) = V^iqxW^iqy) - (y^(q))^ 



(3.4) 



and Ky{q) = KxiPq), Vy{q) = v^iPq) where Pq = 
P{qx,qy) = {qy,qx)- From Eq. (3.3) for G~\ we can 
calculate 



G 



where 



V = 



1 / ^[^ + -y^l 



V 



, (3.5) 



+ (3.6) 



is the determinant of G~^. 

In order to determine the dimensions of operators, we 
calculate the leading dependence of correlation functions 
such (^^(r, f)) on system size L. Thus we can consider 
the function 

//,2. ,^^ _ I dq^dqydw + %gy) 
l^^x(r,t))-y -^2^^^ ^ (3.7) 
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The leading L dependence is related to the infrared di- 
vergence of the integral just introduced. This infrared 
divergence comes purely from the integration over qx.oj. 
We can write the integrand as Kx/{^ — I" pl^s a re- 

maining part, the integral over whicfi is free of infrared 
singularities (see Appendix C for details). Thus, it is 
easy to see that the leading L dependence goes as 



dqy 



(3.8) 



where K{qy) = Kx{0,qy). Notice that this is precisely 
what we had for a single array of parallel wires. A sim- 
ilar analysis yields {9y) = (6*^). Also note that cross- 
correlations of the form 9x-0y are finite as L goes to in- 
finity. 

We also need to consider correlation functions in the 6 



variables. To do so, we start with the action of Eq. 3.1 



and integrate out the 9 variables. The effective action as 
a functional of the 6 variables is 



Sa, — — 



duodqxdqy 
2 / 2^ 



with a = x,y; b = x,y. Here 
V'qxqy 



V^qxqy 



(3.9) 



(3.10) 



where 



(3.11) 



Note that k is different from n defined for the 9 correla- 
tion functions. From G7^, we calculate 



G„ 



where 





-V^qxqy 



-V^qxqy 



Vxql 



(3.12) 



-{V^lWy 
1-1 



(3.13) 



is the determinant of G|^ . G^ can be used to calculate 
dimensions of operators involving 0. For example the 
expectation value 



(0l(r,t)) 



dqxdqydoj 



-(* 



Vyql 



(27r)3 



V 



(3.14) 



where again we take ^ < \qx\, \qy\, In the integral, 
as before, the infrared divergence comes purely from the 
integral over qx,0J. Once more, the infrared divergent 
part goes as 



{4>l)^irHL) 



dqy 



27r K,{q; 



(3.15) 



V! 



where n{qy) = K.x{Q,qy) is the same function appearing 
in {91), Eq. (3.8). 

Thus, correlation functions for 9x and 9y can be calcu- 



lated directly from Eq. (^j) . 9x-9y cross-correlations are 
non-singular, whereas, 6x-9x and 9y-9y correlations have 
singular parts with exactly the same functional forms as 
they have in the absence of coupling between layers, but 
with the ti{q) function in expressions for the scaling ex- 
ponents replaced by 



n{qi_) = Hx{Q,qi_) Ky{q^,Q). 



(3.16) 



The same holds for 0-0 correlation functions. Thus corre- 
lation functions within a given array have the same func- 
tional form as for — Q but with different definitions 
of K. Other than renormalizing «((/), the coupling „ 
between the two arrays leaves the dimensions of all oper- 
ators unchanged. This means that it is possible to choose 
interchain interactions within the X and K-grids so that 
these grids form 2D anisotropic sliding phases even in 
t he prese nce o f the inter-grid coupling V^'^ „. Equations 
( |3.2[ ) and ( 3.16| ) define a 2D non- Fermi liquid with scaling 
properties to be discussed in the next section. 

First, however, we must verify that it is possible to 
choose potentials so that this 2D non-Fcrmi liquid is 
stable with respect to perturbations. All pairwise cou- 
plings within a given array, i.e. S'gc,™, <S'cDW,n' ^Icn and 
S'cDw n defined as obvious generahzations of Eqs. (2.19), 
can be rendered irrelevant by choosing K{qi_) as in the 
case of an individual array. We must also consider 
Josephson and CDW couplings between the two arrays, 
which operate at the points of crossing {x,y) = [na, ma) 
of wire m in the X-array and wire n of the K-array re- 
spectively. These take the form 

^SC ^Y. d'^'^^^ cos[y/2^{9x.^m{na) - Oy^nima))] 
ScDW ^ X! y 'i'^'^^^ cos[V2TT{(t)x,m{na) - 0;,,„(ma)) 



-\-2kp{ma — na)]. 



(3.17) 



The dimensions of the cosine operators in the integrands 
are, respectively, 

dq 



Ascoo ^ J ^n{q) = K 
r'^ dq 1 



iCDW,oo 



1 1 



-TT 27r K{q) 



(3.18) 



where we assume that K{q) has the form given by Eq. 
(3.17), A is defined as before, and C = K"{ko)/2K. If k 
is chosen such that Eq. ( |3.16 ) is satisfied for each array, 
then S^q' and S^^^ are automatically irrelevant. Thus, 
we do not need any further fine tuning of k to get a stable 
CSLL phase. 
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Having established a region of stability of the CSLL 
phase, we now investigate the nature of the correlation 
functions. Consider once more the correlation function 

(3.19) 

which corresponds to the component of the density- 
density correlation function modulated at 2kp. In the 
absence of terms such as Sq^^j^j, this correlation function 
vanishes for y 0. Thus, though irrelevant, the presence 
of changes the nature of the correlation functions. 

In its presence, to lowest order in V''^^ , we obtain 



J2k]p{x+y) 



dxi [ 



l^^iV2TT(4>^{x,y=ma)-tp^{xi,y=ma))j^ 

where ()o is the expectation value with respect to the 
CSLL fixed point. We need the asymptotic form of the 
correlation function for large x, y. Following Eq. (2.12) 
we obtain, for example. 



^giV27r(</,,(x,O)-0,(O,O))^^ 



rAcE 



Similarly 



/giV27r(0j(x,y)-0x(x,O))\ 



(3.21) 



(3.22) 



Note that in the slidin g ph ase, AcDW.n is greater than 
unity and thus in Eq. ( ^.20 ), the largest contribution to 
the integral comes from xi close to and to x. It thus 
follows that for large x and y, goes as 



G 



X C cos[2kF{x + y)] 



{xy) 



Ace 



(3.23) 



where C is a constant that is proportional to (V"^^)^. 
Thus the correlation function decays as a power law 
in all directions. Notice that G^ is not isotropic but 
exhibits a square C^y symmetry. 

The stability of the CSLL phase for the gapless case 
follows along the same lines. If there are no marginal 
inter-array spin-dependent coupling terms, then k„ ~ \, 
and we define a renormalized Kp{q±). The stability of 
the CSLL phase is identical to the stability of the smec- 
tic phase on a single array, with a fixed-point action de- 
scribed by the renor mali zed f uncti on Kp(q±). Also, pro- 
ceeding as in Eqs. ( 3.2C ) to ( 3.23 ), we now expect the 
single-electron correlation functions to exhibit power-law 
decay in all directions. 



B. Crossed 2:1 array 

The above analysis can also be extended quite easily 
to a three-dimensional stack of alternate 2D X- and Y- 
arrays. We could also think of such a stack as a three- 
dimension array of wires running along the X-axis, in- 
termeshed with a 3D array of wires running in the in- 
direction. Thus the fixed point action would be of the 
form Sx + Sy + Sxy where Sx and Sy are the actions for 
the 3D arrays formed by wires running along the X-axis 
and y-axis respectively, Sxy represents the inter-array 
Coulomb interactions. Thus the fixed-point action is 

2 J (27r)4 [ly^lzjqxWxl 

+V'{qx,qz)ql\Oy\' 
+VHQy,<lz)<ll\<t>x\' + V^iq^,q,)kl\(l)y\' 

+{V^^ {qx, qy, qz)qxqy(t)x<l)*y + c.c.} 

-iujq^ieicj)^ + c.c.} - iuqy{9*y(j)y + c.c.}] (3.24) 

where (j)x,y and 9x,y are functions of qx, qy, and qz, 
and {qx,qy,qz) represents the interactions between 
the X- and K-arrays. For interactions only between 
nearest neighbor layers, we obtain {qx,qy,qz) = 

V'^{qx,qy){l + e*^^). Integrating out the (j) variables, we 
are left with an effective action 



„ I [ dujdqxdqy 
^'~2j^2^ 



1 



Vxi(l)ql 



Vxiq) 



1 



Ky{q) 

rXY 



Vy{Ci) 



{F/^(q)^2M;+c.C.}], 



\ex 



(3.25) 



where 



Kx(q) 



7(q) 



V<'{qy,qz)V^{qx.qzy 
V«(g„g.)7(q) 



V-^iqx^qz) 



7(q) 



l{ci)^V*{qx)V^{qy)~\V''^\^ 



(3.26) 



and Hiy{q) = Kj;(Pq), Vy{q) = Vx{Pq) where Pq = 
P{qx,qy,qz) = {qy,qx,qz)- Proceeding exactly as in the 
previous case, we find that 



q2\ _i rT\ I '^ly'^lz 
(2^)2 



7rln(L) 



Kx{0,qy,qz) (3.27) 



and 



dqydqz 1 
(27r)2 Kx{Q,qy,qz)' 



(3.28) 
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The stability of the three-dimensional crossed stack is 
precisely the same as the stability of a three-dimensional 
stack of parallel quantum wires with the Luttinger Liq- 
uid parameter K{qy,qz) set equal to Kx{0,qy,qz) of the 
crossed stack. As before, there are no additional singu- 
larities due to the coupling between the crossed arrays. 



IV. TRANSPORT PROPERTIES 

We now investigate the transport properties of the slid- 
ing Luttinger liquid phases. The conductivities of an ar- 
ray of parallel wires has been considered by Emery et. 
al In a pure system the conductivity along a wire 
is infinite. In the presence of impurities, the resistivity 
along the wires vanishes as pa]: 



with 



iCDW,oo 



(4.1) 



(4.2) 



The conductivity perpendicular to the wires for an array 
of parallel wire can be calculated ||l^,|2^ using the Kubo 
formula, giving us 



(4.3) 



with (Ref. |£8[) a±^ — 2Agc — 3, where Asc is the min- 
imum of Age 1 and Age, 2 (For details see Appendix D). 
The conductance, Cc, arising from the Josephson cou- 
pling at the contact between the crossed wires, can be 
calculated similarly using the Kubo formula, and satis- 
fies 



(4.4) 



where ac = 2Agc — 3. In this section we focus on 
the gapped case. In the gapless case, p\\ , a± and Cc 
still exhibit power-law behavior even though the major 
contribution to perpendicular conductivities may come 
from single-particle hopping. 

Thus we can model our 2D non-Fermi liquid as the re- 
sistor network depicted in Fig. (7) with nodes at the ver- 
tical Josephson junctions between the arrays at (x, y) — 
{na, ma). The nodes of the A'(y)-array are connected by 
nearest neighbor resistors with conductances (T|[ ~ pj'^ 

if they are parallel to the x(y)-axis and a± if they are 
perpendicular to the x-axis(2/-axis). Nearest neighbor 
nodes of the X and F-arrays are connected by resistors 
of conductance ac- In the continuum limit, the 2D cur- 
rent densities in the plane of the a grids (a = X, Y) is 



af^Ef where 
''J J 



X 



(T|| 

cr_L 

(T_L 

an 



(4.5) 



(4.6) 




FIG. 7. A schematic depiction of the 2D non-Fermi liquid 
as a resistor network, with two parallel arrays of wire running 
along the x and y-axes, with nodes in the z direction 



and E" is the in-plane electric field in plane a. The 
current per unit area passing between the planes is J„ = 
{ad a^){V-^ — V^) where V is the local voltage. In this 
limit, the local voltages satisfy 



X 



(4.7) 



where and are current densities (current/area) 
injected, respectively, into the X and F-grids. If no 
currents are injected, then this equation is solved by 
= = — E • X to produce a total in-planar cur- 
rent density 



Ji — J'i 



)E^ = {a«+a^)E,. (4.8) 



Thus under a uniform electric field, the double layer be- 
haves like an isotropic 2D material with in-plane conduc- 
tivity a = a\\-\-a^_ ~ a\\ , or equivalently with an isotropic 
resistivity that vanishes as py ~ T"ii . 

We could also consider currents that are spatially 
nonuniform, as they are, for example, when current is in- 
serted at one point and extracted from another. In that 
case, there is a crossover from isotropic to anisotropic 
behavior at length scale 




rp-{a\\+a^)/2 



(4.9) 



that diverges as T ^ 0. To illustrate this crossover, we 
calculate explicitly the case where a current / is inserted 
at a point ri on the X-array and extracted at another 
point r2 on the X-array. Then 



=/[<5(r-ri)-<5(r-r2)]. 



(4.10) 



Using equations (4.7) and (4.10), one can solve for the 



resistance between these two points: 
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d^q 1 



oiq(r2-r2) 



(27r)2 



.9(q) 



(4.11) 



where 



5(q) 



f*[(^ii 



ml 



For 



^ -<j±ql-(^\\ql 



(4.12) 



(4.13) 



(7(q) takes the simple form (o-|| + cr_L)g^. If |ri —V2\ ^ I, 
with Z de fined in Eq. ( [4.9| ) , then the integral over g in 
Eg. (|4.11 ) is dominated by small q satisfying equation 
(4.13). Thus for |ri — » /, the system has approxi- 



mately the same resistance as an isotropic conductor with 
conductivity (T|| . If we inserted current at a point on 
the X-array and ext racte d it from the y-array we would 
have the form Eq. ( 4.11 ) for the resistance with a dif- 
ferent function q (q) whose small q limit is still given 
by (ay -I- a-±)q^. Thus for |ri — ^ / the resistance 
is approximately independent of whether the current is 
inserted into (or extracted from) the X or the F-array. 
More generally, in a region where T = 0, by inspection of 
equation (4.7) it can be seen that inhomogeneities in the 
voltage (or difference between V'^ an d V ^) would heal 
over the lengthscale I as defined in Eq. (|4.9| ), and at longer 
lengthscales the system would behave isotropically. This 
length diverges as T ^ and at T = 0, current can only 
be carried along the wires: the resistance between wires 
in a grid or between grids is infinite. 

We could, in addition, investigate the frequency- 
dependent zero-temperature conductivity. By arguments 
similar to those leading Eqs. (4.1) to (4.4), we obtain 



(4.14) 



Q!j[, a±, and ac are the same as before, though the coef- 
ficients are now different (and complex, in general). At 
finite u, the long- wavelength resistivity is isotropic as be- 
fore, and vanishes as (w) ^ . 

We could also consider extensions of these calculations 
to three-dimensional stacks of crossed arrays. As we saw 
in the previous section, it is possible to get a stable sliding 
phase in such a system. The conductivity now has a 
three-dimensional character, with conductivity along the 
planes given by (J\\/d, d being the separation between 
adjacent X-arrays, but conductivity in the third direction 
given by ad a. Thus the conductivity along the planes is 
much larger than the perpendicular conductivity. 



V. CONCLUSION 

In conclusion, we have demonstrated the existence of 
non- Fermi metallic phases in two and three dimensions, 
that are stable with respect to a wide class of perturba- 
tions. We consider both spin-gapped systems and gapless 
systems which exhibit spin-charge separation. Our cen- 
tral results pertain to the stability and properties of the 
crossed sliding Luttinger Liquid (CSLL) phase. This is 
a remarkable phase, which could be identified as a two- 
dimensional Luttinger liquid. The correlation functions 
in this phase exhibit power-law decay along the planes, 
and the finite-temperature long-wavelength electric con- 
ductivity, which is isotropic along the planes, diverges as 
a power-law in temperature T as T — > 0. The importance 
of this work is that it provides a perturbative access to 
non-Fermi liquid fixed points in two and three dimen- 
sional systems, something that has proven to be quite 
difficult in the past . This work could be of significant 
relevance for higher dimensional strongly correlated elec- 
tron systems in general, and to the normal conducting 
phases of the cuprates in particular. 

RM and TCL acknowledge support from the National 
Science Foundation under grant DMROO-96531. We ac- 
knowledge useful discussions with S. Kivelson, E. Fradkin 
and A. Vishwanath. 



APPENDIX A: 

In this appendix we sketch out the steps leading to the 
asymptotic form for the correlation function G^{x, r = 0) 
for large x (see Eqs. (2.12) and (2.13)). 



= e-2-(W.(-)-*.(0)f > cosi2kpx). 



c.c. 



(Al) 



It can be easily checked that 

((</>,(x)-0,(O))2)=2 
1 — cos(a:;g|| ) 



dq^\dq±duj 
(27r)3 



«(9±)(t'(9±)'zf + ^jf^) 



(A2) 



Next we carry out the integration over gy and uj obtain- 
ing, for large x, 



dqnduj 



1 — cos(a;(7|| ) 



7rlog(a::) 

^{q±){v{q^)ql + ^) ^ <q±) 



Fiqj 



(A3) 



where F{qj_) is some function of q± which depends on k, 
v{q±) and the momentum cutoff. From Eqs. (A.l) and 
(A. 3) it follows that 
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((0,(x)-0,(O)nc 

Using this, we obtain 

G4,{x,Q) 

where 

AcDW.oo 



dq±_ 1 



Ai coa{2kpx) 
dq± 1 



log(a;) + const.. 

(A4) 

(A5) 



Eqs. (2.14) and (2.15) follow along similar lines. 



(A6) 



APPENDIX B: 

Here we outline how the integral 

r dq (1 -cosng) 
^ i-. ^ /(.) ' ^""'^ 

needed to calculate AcDW,nj can be solved exactly. Here 
/(g) = l + Aicos(g) + A2Cos(2(7), (B2) 



using expression ( 2.30 ) for Ai and A2. We can rewrite 



/((?) = 2 A2 (cos g + u+) (cos g + u ) 



where 



and 



2A2 



±iD 



A? 



2A 



n-' — L. J 2 — 

4A2 + A2 ^ " A2 ■ 



It is easy to check that 



In = - 



V2AA2 



zlmJj 



(B3) 



(B4) 



(B5) 



(B6) 



with 



7+ 

J2,n 



dq 1 — cos{nq) 
27r M+ + cos q 



— :[Jl,n + J2,n] 

m 

dz 

1 + 2u+z + z2 
I f dz [1 + 



2 / z" 1 + 2u+z + z2 



(B7) 



where z = e"', and the z integral is over the unit circle 
centered about the origin. The integrands have poles at 
z = {0,z^,Zq): 



Using (A8), it is easy to check that Zq 



(B8) 



Thus, either Zq or Zq lies inside the contour of integra- 
tion (the unit circle). Using the method of residues, it is 
now straightforward to calculate the integrals. We sim- 
ply need to sum over the residues of the poles enclosed 
within the contour of integration. In order to express our 
results, we distinguish two cases: 

1. \z+\ < 1. Then 



J. 



+ _ 



(^0^ - ^0 ) 

n-l 



1 - 



(i + (4r) 

2(4)" 



E 



1 



- (z+)™+i(zo-)- 

^ -[1-iz^r] 



(4 - ^0 ) 



(B9) 



2. Izn I < 1. In this case 



/+ - 



E 



2(^0")" 



'„(z+)™+i(zo-)«- 

[i-(4n- 



APPENDIX C: 



(BIO) 



Here we consider the infrared divergence of the integral 



/ = 



dqxdqydu! \vy ~^ "^vly 
(27r)3 1^ 



(CI) 



where 



V = 



— + — + ^^9, - {V^y^- (C2) 



At first sight, it may be appear that the integral is diver- 
gence free, since by power counting, there are two powers 
of Q (where Q = (uj,qx,qy)) in the numerator multiply- 
ing d^Q, and four powers of Q in the denominator. This 
seems to indicate that the integral is finite as L — > 00. 
Notice, however, that if is set equal to zero in V, the 
integral can be written as 



I 



dq^dqydw K^(g 



27r3 



Vxqi 



(C3) 



which is clearly infrared divergent. This divergence 
comes purely from the integration over qx,^- It turns 
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out that even in the presence on the divergence conies 
purely from the integration over q^^ to. To obtain the in- 
frared divergent part we write 



Kn-V 



R 



(C4) 



where R is the remaining piece, and our task is to show 
that its integral has no infrared divergence. Let us write 
R = Ri + i?2 where 



i?2 



The step function Q{t) may be written as (l-|-Sign(i))/2 
where Sign(i) is +1 for positive t, and —1 for negative t. 
In the spin-gapped case the contribution to the param- 
agnetic part comes from superconducting pair-hopping. 
The paramagnetic and diamagnetic terms can be com- 
bined to give 



o-_l(cj) = — 



dx / dte'^^^l + Sign{t)) 



x[n>(x,<) -n<(a;,i)] - = 0)] (D3) 



where 



n> = -i 



(J(a;,j,i)J(0,0,0)) 
^(J(0,0,0)J(a;,j,t)). 



(D4) 



1 

V 



Vxq: 



(C5) 



The integral of Ri has no infrared divergence. To check 
that this is true for R2 as well, we note that in the ex- 
pression for i?2 the term in the square brackets can be 
written as 



f^x f^y ^^/^ 



(C6) 



Now, by noticing the powers of it is easy to see 

that the integration of R2 has no divergence. Thus the 
infrared divergent part of / can be written as 



dqxdqydu Kx{0,qy) 



(27r)3 ^ 
Equation (3.8) now follows easily. 



Vxql 



(C7) 



APPENDIX D: 

We demonstrate explicitly the leading dependence of 
the perpendicular conductivity on temperature. Accord- 
ing to the Kubo formula the transverse conductivity is 
given by 



(J±{uj) = — 



n±(a.) 



(Dl) 



Since II^ — II^ is odd in t, the real part of the conduc- 
tivity is given by 



crj_(tj) = - 
U! 



dx / dte"^'[U>{x,t)^U<{x,t)] 



(D5) 



Note that the DC transverse conductivity is purely real, 
and can be obtained from i7j_ (lo) by taking the limit 
uj 0. is related by analytic continuation to the 
Matsubara correlator IIm{x,t) in the upper half plane 
of complex t space, and 11^ is related to II]\i{x,t) in 
the l owe r half plane. Thus we could view the integral in 
Eq. ( |D5| ) as an integral over the Keldyish contour shown 
in figure 5a. This contour can be distorted to the contour 
shown in figure (5b). Note that n>{t+§) = n<{t~'-§) 
where /3 — -^Kp . Thus, we obtain 



Cr_L(w) = - 



dx 



dt 



--sinhf'^f 



n>(x,i + ^) 



(D6) 



The next step is to calculate the Matsubara correla- 
tor T\m(x^t). To begin with, we only consider nearest- 
neighbor hoppings between wires. To lowest order in 
and T : 

nM(a;,r) =^(Ji(a;, J, T)Ji (0,0,0)), (D7) 



where 



where the first term represents the paramagnetic contri- 
bution with 



n±(t^) 



dx 



dm{t)e' 



([Ji(a;,j,t),Ji(0,0,0)]) 



(D2) 



being the retarded current-current correlator, and the 
second term represents the diamagnetic contribution. 



aJi[s\ii{9j{x, t) 
+ sin(6'j+i(x,T) 



Vi(a;,r)) 
9,(x,t))]. 



(D8) 



Here a is the distance between adjacent wires, and the ex- 
pectation value is taken with respect to the sliding fixed- 
point Lagrangian. The correlator can be written as: 



= J2a2cxp[-/(a;,r)] (D9) 
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a) 



p/2 



where 77 = Asc,i- The a; — > hmit of the above expres- 
sion yields 



a;(a; = 0)~r(2^^'=.i-3). 



(D14) 



When the velocity w is a function of q±, the integral 
can no longer be solved exactly. However, the leading T 
dependence of the conductivity remains unchanged. To 
check this, we follow the previous set of steps and arrive 
at the expression 



dxdt 



n 



:t;(gi)cosh(f +t)cosh(f -i)]^''^ 



(D15) 



-p/2 

b) 

FIG. 8. a) The Keldyish contour in complex t-plane, with 
real time along the X-axis. In b) we depict how the contour 
is deformed in order to evaluate the integral. 

where 



(diix, t) - e, (0, o)){e„{x, r) - Ooio, o))) 

/ , I dq±dqx[{0{qx,q±,u})0*{qx,q±,oj)} 



x(l - cos{qxX + ljt)) X (1 - cosgj,)]. (DIO) 

Let us first consider a simpler case, where the velocity 
Vs has no dependence with q±. Then 



ILM{x,T) = a^J^ 



{TTTaJvf^ 



[sinh(7rT(a;/t; + ir)) sinh(7rT(a;/i; - ir))]'' 



(Dll) 



where 77 = (5sc,i) and is the spatial cutoff along x. 
Thus we may write 



I 



sinh ( ^) 

3h(.T + i) cosh(£ - t)Yi 



X / dxdt- 



(D12) 



where x = vxT jv and t = irtT. By introducing new 
variables x+i and ^ — t, we carry out the above integrals, 
giving 



2 r„n„^2r, !^^^|i(i^^„^^T.^2^-3 

8uj/T 



-p2 (V I \ -r2 (n 



(D13) 



where x = wxT, t = wtT and J2qj_ '^g± = V = -^SC,!- 
The result can be expressed in the scaling form 



(D16) 



In the limit oj ^ 0, the integral is finite and T indepen- 
dent, implying F{()) is finite. Thus ct_l(w = 0) ~ T", 
where rj ~ Asc.i for nearest-neighbor hopping. In ex- 
actly the same manner, we can calculate the contribu- 
tion to aj_ from next-nearest neighbor pair-hopping. It 
has the same scaling form as before with Asc,i replaced 
by Asc,2- 
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